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ABSTRACT 

Bimodal globular cluster metallicity distributions have now been seen 
in a handful of large ellipticals. Here we report the discovery of a bimodal 
distribution in the dominant group elliptical NGC 5846, using the Hubble 
Space Telescope's Wide Field and Planetary Camera 2 (WFPC2). The 
two peaks are located at V-I = 0.96 and 1.17, which roughly correspond 
to metallicities of [Fe/H] = -1.2 and -0.2 respectively. The luminosity 
functions of the blue and red subpopulations appear to be the same. 



^Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555 
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indicating that luminosity does not correlate with metallicity within an 
individual galaxy's globular cluster system. 

Our WFPC2 data cover three pointings allowing us to examine the 
spatial distribution of globular clusters out to 30 kpc (or 2.5 galaxy 
effective radii). We find a power law surface density with a very fiat 
slope, and a tendency for globular clusters to align close to the galaxy 
minor axis. An extrapolation of the surface density profile, out to 50 
kpc, gives a specific frequency S^v = 4.3 ±1.1. Thus NGC 5846 has a 
much lower specific frequency than other dominant ellipticals in clusters 
but is similar to those in groups. The central galaxy regions reveal some 
filamentary dust features, presumably from a past merger or accretion of 
a gas-rich galaxy. This dust reaches to the very nucleus and so provides 
an obvious source of fuel for the radio core. We have searched for proto- 
globular clusters that may have resulted from the merger/accretion and 
find none. Finally, we briefiy discuss the implications of our results for 
globular cluster formation mechanisms. 
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1. Introduction 

One of the most interesting recent develop- 
ments in globular cluster (GC) studies is the 
discovery that at least some elliptical galax- 
ies have a bimodal GC metallicity distribu- 
tion. A bimodal distribution is also seen in 
the Milky Way GG system. These bimodal, 
and sometimes possibly trimodal, distribu- 
tions are more metal-rich in the mean than 
those of the Milky Way but the same general 
conclusion can be drawn from their existence; 
galaxies with multimodal metallicity distribu- 
tions have not undergone a simple monolithic 
collapse but have experienced episodic star 
formation. To date, the best examples of mul- 
timodal GG metallicity distributions are cD 
galaxies with extremely rich GG systems (i.e. 
with specific frequency S^v ~ 15) such as M87 
in Virgo (McLaughlin, Harris & Hanes 1994), 
NGG 1399 in Fornax (Bridges, Hanes & Har- 
ris 1991) and NGG 3311 in Hydra (McLaugh- 
lin et al. 1995). Two galaxies with 'nor- 
mal' specific frequencies have also revealed bi- 
modal distributions i.e., NGG 4472 with Sat 
= 5.6 (Geisler et al. 1996) and NGG 3923 
with Sat = 6.4 (Zepf, Ashman & Geisler 1995). 
Thus the formation mechanism responsible 
for two distinct GG populations must be able 
to explain both the populous GG systems in 
cD galaxies, and the 'normal' richness GG sys- 
tems in ellipticals. 

An interesting galaxy in this regard is NGG 
5846. Although it dominates a small, com- 
pact group of galaxies it has an average or 
slightly below average specific frequency. De- 
tection of a bimodal GG color distribution in 
this galaxy would add a further constraint to 
GG formation mechanisms. The galaxy is not 
a cD but a giant EO elliptical, with an effective 
radius of 11.6 kpc (Bender, Burstein & Faber 
1992) and an absolute magnitude of My = 



-22.6 (Faber et al. 1989), assuming a dis- 
tance modulus (m-M) = 32.3 from Forbes, 
Brodie & Huchra (1996; hereafter Paper 1). 
This would place NGG 5846 some 13 Mpc 
more distant than the Virgo cluster. Harris & 
van den Bergh (1981) estimated the total GG 
count, from a photographic study, to be 2200 
± 1300. For My = -22.6 this corresponds to 
Sat = 2.0 ± 1.2. (Note that Harris & van den 
Bergh used My = -21.7 and hence derived Sat 
= 4.5 ± 2.6.) 

The inner region of NGG 5846 contains 
some dust and a compact radio core with 
a possible extended component (Mollcnhoff, 
Hummel & Bender 1992). The X-ray emis- 
sion from hot gas extends ~ 100 kpc from 
the center of the galaxy, encompassing two 
other early-type galaxies in the group (Bier- 
mann, Kronbcrg & Schumutzlcr 1989). Bier- 
mann et al. calculated the hot gas mass to 
be ~ lO^^M©, and the total virial mass to be 
almost 1O^^M0. Recent ROSAT mapping re- 
veals that the soft X-ray emission is spatially 
coincident with the Ha emission and dust 
morphology (Goudfrooij & Trinchieri 1996). 
The galaxy shows many characteristics simi- 
lar to those of M87, which lies at the center of 
the Virgo cluster. Both arc luminous galaxies 
(My = -22.5 for M87) which dominate the 
local potential and lie at the center of an ex- 
tensive X-ray envelope. The virial mass and 
coohng flow infall rate are also similar for the 
two galaxies. The central velocity dispersion 
is 361 km s^^ for M87 and 278 km for 
NGG 5846 (Faber et al. 1989). They differ 
greatly in one important regard, namely GG 
specific frequency (Sat = 14 ± 0.5 for M87). It 
has been speculated that some of M87's GGs 
were formed as a result of its location at the 
center of a large gravitational potential, ei- 
ther from the coohng flow (Fabian, Nulsen & 
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Canizares 1984), or by accretion from neigh- 
boring galaxies (e.g. Forte et al. 1982) or 
from the field (West et al. 1996). Knowl- 
edge of the GC system properties of NGC 
5846, particularly in comparison with giant 
ellipticals in different environments, may pro- 
vide important clues about globular cluster 
and galaxy formation histories. 

In paper I we presented the V band WFPC2 
data for NGC 5846. The GC luminosity func- 
tion (GCLF) was found to be the same in each 
of three separate pointings, ranging from ~ 
3 to 30 kpc from the galaxy center. While 
taking into account incompleteness, photo- 
metric errors and background contamination 
we fit the combined GCLF using the maxi- 
mum likelihood technique of Seeker & Harris 
(1993). This gave a turnover magnitude of 
my = 25.05 ± 0.10 and a dispersion of cr = 
1.34 ± 0.06. Here we include the WFPC2 I 
band data and discuss the color and spatial 
distributions of the GC system in NGC 5846. 

2. Observations and Data Reduction 

Details of the observations and data re- 
duction are described in Paper I. Briefly, we 
obtained WFPC2 images of NGC 5846 at 
three pointings, in two filters. The point- 
ings are located ~ 1.5 arcmin north of the 
nucleus, ~ 2.5 arcmin south and one with 
the galaxy centered in the PC CCD. For each 
filter, we have a pair of images giving a to- 
tal of 2400 seconds in F555W and 2300 sec- 
onds in F814W. Objects were detected us- 
ing DAOPHOT (Stetson 1987) with a S/N 
threshold of 4 per pixel. Aperture magnitudes 
(in a 2 pixel radius) were converted to total 
Johnson V and Cousins I mags. A Galac- 
tic extinction correction (Ay = 0.11 and A/ 
= 0.05) was also apphed. The DAOPHOT- 
selected object list was checked for matches to 



hotpixels and all sources with FWHM greater 
than 3 pixels were removed. Contamination 
in the final object list is expected to be less 
than 5%. Completeness tests for the V band 
data indicate a 50% incompleteness at V = 
25.95. Photometric errors at this magnitude 
are 0.16. In this paper we include GCs in 
the central PC which are located at galacto- 
centric distances of greater than 7 (1 kpc) 
i.e., beyond the strongest dust lanes. Globu- 
lar clusters located interior to this radius are 
discussed separately in section 4.1. 

3. Results 

3.1. The Galaxy 

We have subtracted a smooth model of the 
galaxy from the PC CCD image of the cen- 
tral pointing. A grey scale of the residuals 
is shown in Fig. 1 (Plate **) which reveals 
several GCs along with filamentary dust fea- 
tures. The DAOPHOT selected GCs, beyond 
7 of the galaxy center, are indicated. In Fig. 
2 we give radial profiles of some parameters 
from the I band galaxy model (which is less ef- 
fected by dust than the V band) . The surface 
brightness profiles show that the V-I color 
is roughly constant between 0.1 and 3 kpc. 
Both V and I profiles show a change in slope 
at about 0.3 kpc or 2.3 . The ellipticity pro- 
file confirms that the galaxy is EO-1 in the 
central regions. The position angle changes 
from 70° to ~ 100° near the nucleus, and 4th 
cosine term (see for example Forbes & Thom- 
son 1992) becomes more boxy which indicates 
extra fight at 45° to the major axis. 

3.2. The Globular Clusters 

The color-magnitude diagram for our full 
sample of 837 GCs, from all three pointings, 
is shown in Fig. 3. The GCs range from V ~ 
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21 to V ~ 26. There is a hint of bimodahty, 

best seen between 24 < V < 25, otherwise 
the average color is relatively constant with 
magnitude. 

The mean color, weighted by photometric 
error, is V-I = 1.15 with a rms scatter of ± 
0.2 ± or in terms of an error on the mean of 
± 0.01. In the subsequent analysis, we have 
confined the sample to have colors within 3a 
of the mean. This excludes 63 objects with 
extreme blue and red colors. These objects 
appear to be distributed randomly through- 
out the three pointings, in particular they 
are not centrally located. The objects with 
the extreme colors tend to be the faintest 
sources, suggesting that photometric errors 
have affected their colors. Histograms for the 
color-selected subsample are shown in Fig. 4. 
Here a color histogram is given for each of 
the pointings and for all three pointings com- 
bined. A clear bimodal color distribution is 
seen in the north and south pointings (and 
is reflected in the combined histogram), with 
peaks located at V-1 = 0.96 and 1.17. The 
central pointing reveals the red peak, with a 
hint of some additional blue GCs. 

Bearing in mind that the conversion of 
color into metallicity is uncertain, we use the 
relation of Couture et al. (1990) to derive 
[Fe/H] from V-1. The rms error in this rela- 
tion (~ 0.25 dex) dominates over the error as- 
sociated with measuring the mean GC color. 
Their relation gives a mean metallicity for the 
GC system of [Fe/H] = -0.3 ± 0.25. In Fig. 
5 we show GC metallicity function of 

galactocentric distance expressed in kpc. A 
weighted fit to [Fe/H] versus log of the dis- 
tance in arcsec gives a hnear slope of -0.13 
± 0.03 over the radial range of 1 to 30 kpc. 
This indicates that the gradient is significant 
at the 4(7 level. 



The surface density distribution of GGs 
around NGC 5846 has been calculated using 
two different methods. The first, applied to 
the central pointing, is to make use of the fact 
that out to some radius the WFPC2 gives al- 
most complete 180° coverage for one hemi- 
sphere of the galaxy. This hemisphere is ± 
90° about the WFPC2 V3 axis. Thus we 
have counted up the GCs in six half-annuli 
from 1 to 13 kpc. The number in each half- 
annulus is doubled to give the full 360° cov- 
erage, and divided by the of the annu- 
lus. For the two offset pointings, the orienta- 
tion of the WFPC2 field-of-view with respect 
to the galaxy is much more complicated. In 
the second method, we simply count up the 
number of GCs per CCD and divide by the 
CCD As the GCs are not uniformly 

distributed across the CCD, we use the me- 
dian galactocentric distance as a weighted dis- 
tance. Surface densities from the two meth- 
ods are given in Fig. 6 which shows a smooth 
transition between the two methods. A fit to 
the surface density distribution using 

p = PoT"" 

gives a = -0.7 ± 0.1 and po = 4.8 ± 0.5 
(assuming negligible background contamina- 
tion). 

We can use the surface density profile to es- 
timate the total number of GCs around NGC 
5846. First we must correct for the faint un- 
detected GCs. From our GCLF fitting in Pa- 
per I we estimate that ~ 25% of GCs were 
not detected. Thus the normalization con- 
stant Po should be increased to 6.0. The un- 
certainty in calculating the total number of 
GCs is dominated by the choice of the outer 
radius. The GC counts of Harris & van den 
Bergh (1981) reached the background level at 
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~ 6' (50 kpc). Integrating the profile out to 
this radius gives a total of 4670 GCs and a 
corresponding specific frequency of = 4.3 
± 1.1. Integrating to 60 kpc would increase 
the number of GCs by 1270. In addition to 
a global Sat, we can also estimate S at as a 
function of radius from our central pointing. 
Using the integrated galaxy V magnitude at 
various radii, from the catalog of Longo & 
de Vaucouleurs (1983), we can calculate the 
'local' Sjv value from our GC counts. The S^v 
values as a function of galactocentric distance 
are shown in Fig. 7. 

The radial distribution of V magnitude is 
shown in Fig. 8. The mean magnitude shows 
little change with distance from the galaxy 
center. The azimuthal distribution of GC 
counts is shown in Fig. 9. Here we have re- 
stricted GCs to lie within 13 kpc and between 
position angles of -25° and 155° in order to 
obtain complete coverage in one hemisphere. 
We find a slight enhancement of GCs around 
20° and a deficit around 100°. 

4. Discussion 

4.1. Central Regions 

The galaxy central regions reveal a num- 
ber of GCs and filamentary dust. The central 
dust indicates a past merger or accretion of 
a gas-rich galaxy and appears to reach right 
into the nucleus. Given the short dynami- 
cal timescale, for galactocentric distances of 
< 0.5 kpc, the dust (and gas) provides an 
obvious source of fuel for the compact radio 
core (MoUenhoff et al. 1992). A number of 
proto-globular clusters have been detected in 
currently merging galaxies (e.g. Holtzman et 
al. 1992; Whitmore et al. 1993; Whitmore 
&; Schweizer 1995; Holtzman et al. 1996), 
as predicted by the Ashman & Zepf (1992) 



merger model. These objects, initially blue 
and very bright, will fade and redden with 
time. Thus depending on the epoch of forma- 
tion of GCs, we may expect to see a popula- 
tion of very blue or very red GCs. A search for 
proto-globular clusters in NGC 5846 is prob- 
lematic due to the dust. As mentioned in sec- 
tion 2, we excluded regions within 7 (1 kpc) 
of the galaxy nucleus from our object list. We 
have decided to visually inspect the inner 7 
in both the V and I images for GCs. Although 
not objective, the eye is good at recognizing 
the brighter sources. We find 9 objects, some 
on the edges of dust lanes. The V magnitudes 
and colors range from 23 < V < 24 and 0.9 
< V-I < 1.7. The mean color is V-I = 1.15 
which is the same as the mean for the full 
sample. Thus we find no evidence for a cen- 
trally concentrated population of GCs with 
extreme colors. Furthermore, if we assume 
that all objects are instrinsically redder than 
V-I = 0.9, then we can place a lower limit on 
their age from stellar population models. Us- 
ing Bruzual & Chariot (1996), a cluster with 
V-I > 0.9 and solar metallicity is at least 1 
Gyr old. This suggests that i/ these central 
GCs formed from a gas-rich merger, then it 
occurred over 1 Gyr ago. The spectroscopic 
study of NGC 5846 by Fisher, Franx & Ilhng- 
worth (1995) would tend to support this idea, 
as they found no evidence for a young stellar 
population and concluded that the stars were 
uniformly old (~ 15 Gyr). 

The galaxy surface brightness profiles (Fig. 
2) do not continue to rise with the same slope 
in the very inner regions i.e., the profiles 
flatten off. These so-called 'shallow cusps' 
have been seen in HST observations of sev- 
eral galaxies (e.g. Forbes et al. 1995; Lauer et 
al. 1995; CaroUo et al. 1996). If we fit the V 
band profile with the dual power-law of Lauer 
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et al. we derive an inner slope of -0.2. For 
My = -22.6, this value lies on the correlation 
of inner slope versus absolute V magnitude 
(Forbes et al. 1995). Although the origin for 
this correlation is still subject to debate, it is 
not the result of dust obscuration but rather 
some fundamental scaling relation in galaxies. 
NGC 5846 obeys this scaling relation. 

4.2. Color Distribution 

The distribution of GC V-I color reveals a 
bimodal structure, with peaks around V-I = 
0.96 and 1.17 (Fig. 4). A KS test confirms 
that the blue subpopulation is statistically 
significant and a reduced goodness-of-fit 
indicates that two Gaussians (with six free 
parameters) are a better fit to the color dis- 
tribution than one Gaussian (with three free 
parameters). Assuming the color-metallicity 
relation of Couture et al. (1990), these colors 
correspond to metallicities of [Fe/H] = -1.2 
and -0.2 respectively. In order to study the 
blue and red subpopulations separately, we 
have defined a blue sample of 232 GCs with 
colors of 0.84 < V-I < 1.06 and a red sam- 
ple of 361 with 1.06 < V-I < 1.28. The ratio 
of the red to blue samples changes from 1.70 
for the central pointing to 1.34 for the north 
and 1.54 for the south pointings. This tends 
to suggest that the red subpopulation is more 
centrally concentrated than the blue popula- 
tion. We measure a small radial metallicity 
gradient, A[Fe/H]/AlogR (arcsec) = -0.13 ± 
0.03. In the case of NGC 4472, a metaUicity 
gradient of -0.4 in A[Fe/II]/AlogR (arcsec) 
was seen by Geisler et al. (1996) using Wash- 
ington photometry. They concluded that the 
gradient was probably due to the changing 
proportions of the two subpopulations rather 
than a smoothly varying radial gradient. 

Next we examine whether or not the blue 



(metal-poor) and red (metal-rich) samples 
have different luminosity functions. The de- 
pendence of GC luminosity on metallicity 
pertains to various aspects of GC formation 
and evolution i.e. self-enrichment, the ex- 
tent to which massive GCs can retain heavy 
elements from the first generation of stars 
and the ability of proto-globular clusters to 
sweep up enriched material during passage 
near the galaxy center. There is, at present, 
no generally accepted model of GC forma- 
tion and no firm physical basis for assum- 
ing that the GCLF is universal. It is gen- 
erally assumed that the mass-to-light ratio 
of all GCs is the same and, therefore, that 
the mass and luminosity functions have the 
same form. For Milky Way and M31 GCs 
there is no luminosity-metallicity dependence 
(Huchra, Brodie & Kent 1991), but this ques- 
tion remains largely unaddressed for galaxies 
beyond the Local Group. Bimodal color, and 
by inference metallicity, distributions have 
been reported in a number of galaxies (see the 
list of Kissler-Patig 1996) but a comparison 
of the luminosity functions of the two sub- 
populations has so far been attempted only 
for M87. Whitmore et al. (1995) and El- 
son & Santiago (1996) each find evidence for 
bimodal color distributions in their WFPC2 
images of the M87 GC system. They also 
claim that the red subpopulation is on aver- 
age fainter than the blue subpopulation by 
0.13 (Whitmore et al. ) and 0.3 (Elson & 
Santiago). 

In Fig. 10 we show the V band GCLF 
for the blue and red subpopulations of NGC 
5846 compared to a Gaussian fit to the full 
sample in Paper I. Although there is consid- 
erable scatter in each magnitude bin, both 
the blue and red subpopulations appear to be 
consistent with the same Gaussian (i.e. fixed 
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turnover magnitude and dispersion). This is 

confirmed by a K-S test which indicates that 
the two subpopulations are not statistically 
different. In order to test our ability to detect 
an offset in the two GCLFs we have progres- 
sively faded the red sample. These tests indi- 
cate that the two samples become statistically 
different only for offsets greater than 0.3 mag, 
i.e. we can not confirm offsets of less than 
0.3 magnitudes. Elson & Santiago (1996) do 
not quote the statistical significance of their 
0.3 magnitude offset. For NGC 5846, we con- 
clude that the blue and red GC populations 
have similar luminosity functions. There is no 
strong evidence for a luminosity-metallicity 
relation for GCs within a single galaxy's GC 
system. This does not, however, rule out the 
possibility of a GC luminosity-metallicity re- 
lation between galaxies. Indeed such a rela- 
tionship may allow us to understand the con- 
nection between GCLF turnover magnitude, 
GC metalhcity (Ashman, Conti k Zcpf 1995) 
and parent galaxy total luminosity (Seeker & 
Harris 1993). 

4.3. Spatial Distribution 

The radial variations of GC properties may 
constrain theories of GC formation and de- 
struction. In Paper I we showed that the 
GCLF did not vary with galactoccntric ra- 
dius. Excluding the PC CCDs (which have 
a slightly brighter detection limit), we now 
show in Fig. 8 that the mean GC magni- 
tude is roughly constant, out to 30 kpc. If 
the mass-to-light ratio is constant, this im- 
plies that the mean mass of GCs does not vary 
with distance from the galaxy center. A sim- 
ilar result was found by Forbes et al. (1996) 
for a sample of 14 nearby ellipticals. Destruc- 
tion processes, such as dynamical friction, ap- 
pear to be insignificant over the radial range 



probed by our data. This result, and Forbes 
et al. (1996), appear to rule out the model of 
Murali & Weinberg (1996) on GC evolution. 

Fig. 5 shows that the mctallicity (as de- 
rived from V-I colors) has a slight depen- 
dence on galactocentric radius. This gradi- 
ent is generally smaller than that measured 
in other giant elliptical and cD galaxies, but 
these results seem to be consistent with the 
observed trend of metalhcity with specific fre- 
quency (see Forbes, Grillmair & Brodie 1996). 
Gradients are expected in both the Ashman 
& Zepf (1992) merger and dissipational col- 
lapse pictures. However, projection effects 
tend to reduce such gradients (e.g. Seeker 
1996). The Searle & Zinn (1978) picture of 
coalescing subunits, on the other hand, pre- 
dicts no radial mctallicity gradient. 

The surface density of GCs in log space is 
shown in Fig. 6. A fit to all of the data points 
gives a slope of -0.7. If we exclude the inner 
data point, which lies at the radius expected 
for the GC system core size (see Forbes et 
al. 1996), the slope steepens slightly to -0.8. 
Even at -0.8, the surface density slope is, to 
our knowledge, the flattest measured for GC 
systems (Kissler-Patig 1996). We have not 
made any correction for background galaxy 
contamination. We note that Haynes & Gio- 
vanelli (1991) have studied the region around 
the NGC 5846 group. They find both nearby 
and distant galaxies to the north and west of 
NGC 5846. Our data do not cover the west of 
NGC 5846 but do include a northern point- 
ing. Indeed there appears to be an enhance- 
ment of faint galaxies in CCDs 2 and 3 in 
the north pointing. Our selection process, af- 
ter applying the size and color cuts, excludes 
the obvious resolved galaxies. It is possi- 
ble that compact background galaxies remain 
in our object list in the northern pointing. 
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However, the surface density in the north- 
ern pointing is similar to that of the south- 
ern pointing at a given galactocentric radius, 
including that background galaxy contamina- 
tion, if any, is negligibly small. The WFPC2 
geometry makes it difficult to measure the 
galaxy starlight slope, but we estimate it to 
be ~ -1.5. Thus the GC system appears to 
have a much flatter surface density slope and 
is more extended than the galaxy light. This 
situation is common for elliptical galaxies (e.g 
Grillmair et al. 1994). Since the galaxy light 
falls off faster with radius than the density of 
GCs, we would expect the local Sat value (i.e. 
the number of GCs normalized by the galaxy 
light out to that radius) to increase with ra- 
dius. This is indeed seen in Fig. 7, in which 
the local Sjv increases from ~ 0.2 at 1 kpc 
to 2 at 12 kpc. We expect the Siv value to 
increase smoothly to the global value of 4.3 
± 1.1 at a radius of 50 kpc. The total num- 
ber of GCs out to 50 kpc is estimated to be 
4670 ± 1270, which is roughly twice the esti- 
mate of Harris & van den Bergh (1981) from 
photographic plates. We conclude that NGC 
5846 has a much lower specific frequency than 
other dominant ellipticals in rich clusters but 
is similar to dominant ellipticals in groups 
(e.g. Bridges & Hanes 1994). 

The azimuthal distribution of GCs is gen- 
erally consistent with a spherical distribution, 
which is similar to the underlying galaxy light 
i.e., type EO-1. However, there is evidence for 
a slight enhancement of counts around posi- 
tion angle 20° and a deficit at 100°. If this 
effect is real, the GCs appear to be ahgned 
closer to the galaxy's minor axis than its ma- 
jor axis (P. A. ~ 100°). The nearest galaxy to 
NGC 5846 is NGC 5846A, which hes to the 
south (P.A. - 190°). 



4.4. Implications for Globular Cluster 
Formation 

One of main results from this study is that 
NGC 5846 has a bimodal color and, by in- 
ference, metallicity distribution. A handful 
of other galaxies (usually having high lumi- 
nosity and high S^v) also reveal bimodality. 
This indicates that there has been more than 
one epoch of GC formation in these galax- 
ies. A bimodal metallicity distribution is 
a prediction of the Ashman & Zepf (1992) 
merger model. In this model the blue GCs 
are thought to be the metal-poor popula- 
tion from the progenitor galaxies and the red 
GCs are metal-rich ones formed in the merg- 
ing process. The red GCs are expected to 
be more centrally concentrated than the blue 
GCs (which is hinted at in NGC 5846). How- 
ever, this model has some difficulty in explain- 
ing other features. For example, to produce 
a luminosity of My = -22.6 the merging of 
~ 10 L* spirals would be required. It is diffi- 
cult to imagine a clear bimodal GC distribu- 
tion emerging from such a history. The pres- 
ence of filamentary dust in the central regions 
suggests a recent encounter with a gas-rich 
galaxy. The short dynamical times in this re- 
gion (< 10^ yrs) and the fact that the GCs 
appear to be over 1 Gyr old suggest that GCs 
did not form from this event. 

As mentioned in the introduction, M87 and 
NGC 5846 have a similar X-ray luminosity 
and coohng flow rate and yet M87 has about 
three times as many GCs per unit starlight as 
NGC 5846. These results support the view of 
Bridges et al. (1996) that the bulk of GCs did 
not form in a cooling flow. Other mechanisms 
(such as a two-phase collapse or the accretion 
of external GCs) should be considered to ex- 
plain the GC bimodahty in NGC 5846. We 
note that NGC 5846A, a compact E2 galaxy 
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with Mv ~ -18.5, lies ~ 6 kpc south of NGC 
5846. The GCs of NGC 5846A may have 
contributed to those of NGC 5846 via tidal 
stripping and may contribute even more if the 
whole galaxy is eventually accreted. Unfor- 
tunately, NGC 5846A is not covered by our 
WFPC2 images. 

5. Conclusions 

Using deep multi-filter WFPC2 images, 
we discuss the color and spatial distribu- 
tions of the globular cluster system in NGC 
5846. We find a bimodal color distribution 
with peaks at V-I = 0.96 and 1.17 in all 
three WFPC2 pointings. These roughly cor- 
respond to metallicities of [Fe/H] = -1.2 and 
-0.2 respectively. Such a bimodal distribu- 
tion suggests that at least formation episodes 
have contributed to the GC system. The 
red (metal-rich) population may be more cen- 
trally concentrated than the blue (mctal-poor) 
population. The radial metallicity gradient 
is A[Fe/H]/AlogR (arcsec) = -0.13 ± 0.03. 
Both subpopulations have the same lumi- 
nosity functions, within an uncertainty of 
0.3 mag. Thus we find no evidence for a 
luminosity-metallicity relation within an in- 
dividual galaxy's globular cluster system. Pre- 
vious work on the Milky Way and M31 reached 
a similar conclusion (Huchra, Brodie & Kent 
1991). For a constant mass-to-light ratio, the 
lack of a mass-met allicity relation for glob- 
ular clusters, allows several formation mod- 
els to be ruled out (ones in which a metallic- 
ity cooling criterion imprints a characteristic 
mass). 

The surface density of GCs falls off much 
more slowly with radius than the galaxy light. 
The best-fit power law gives a slope of ~ - 
0.75 which appears to be the flattest globular 
cluster systems known. The globular clusters 



appear to be better aligned with the galaxy's 
minor axis than its major axis, although this 
result is of marginal significance. Wc find 
that the specific frequency increases steadily 
with radius in the inner regions and reaches a 
global value of 4.3 ± 1.1 at 50 kpc. Thus NGC 
5846 has a much lower specific frequency than 
other dominant ellipticals in clusters but is 
similar to dominant ellipticals in groups. 

The galaxy isophotes deviate from pure el- 
hpses and become shghtly boxy (extra hght 
at 45° to the major axis) in the central re- 
gions, perhaps due to the filamentary dust. 
The dust is probably the result of a past 
merger or accretion of a gas-rich galaxy and 
provides the fuel for the compact radio core. 
We have searched for proto-globular clusters 
associated with the dust (and gas) and find 
none. 
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Fig. 1.— Grey scale WFPC2 image of the PC CCD (about 30 x 30 ) of the central region 
of NGC 5846. The automatically selected globular clusters beyond 7" (1 kpc) from the galaxy 
center are indicated by black squares. Several radial dust filaments can be seen reaching into the 
nucleus. 
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Fig. 2. — Galaxy properties within 20" (3 kpc). We show surface brightness, eUipticity, position 

angle and the 4th cosine term deviations from a pure ellipse. The I band data are represented by 
filled circles and the V band surface brightness data by open circles. In the central regions, NGC 
5846 has near constant V-I color and eUipticity. The position angle changes from 70° to ~ 100° 
near the nucleus, and 4th cosine term becomes more boxy. 



13 




.5 1 1.5 2 2.5 

V--I color (mag) 

Fig. 3. — Color-magnitude diagram for a sample of 837 globular clusters around NGC 5846. 
Typical error bars are shown on the left. The globular cluster color is roughly constant with 
magnitude. 
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Fig. 4. — Histograms of globular cluster V-I colors (within 3cr of the mean color). The four panels 
show data from the north, south and central pointings, along with the combined sample. A clear 
bimodal color distribution can be seen. 
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Fig. 5. — Radial variation of globular cluster mctallicity. The solid line represents a weighted fit 
to the data. There is a shallow metallicity gradient. The effective radius of the galaxy (re//) is 
also indicated. The typical random error of each data point is ± 0.25. 
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Fig. 6. — Surface density profile of globular clusters. The surface density lias been calculated 
using two different methods (see text for details) which are represented by open and filled circles. 
Error bars represent Poisson statistics. The dashed line is a best fit to all data points, which gives 
a slope of -0.7. 
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Fig. 7. — Radial variation of local specific frequency. The local specific frequency is the number of 
globular clusters normalized by the integrated galaxy light at each radius. The local Sjv increases 
with distance from the galaxy center, reaching a global value of 4.3 at 50 kpc. 
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Fig. 8. — Radial variation of V magnitude. The mean magnitude in 4 bins is indicated along 
with the rms scatter. There is little or no magnitude gradient beyond the central PC CCD (~ 2.5 
kpc), which has a slightly brighter detection limit. 
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Fig. 9. — Histogram of globular cluster position angle. Error bars represent Poisson statistics. 
North and East are indicated. The galaxy major axis lies at ~ 100 . The globular clusters reveal 
a slight enhancement at ~ 20° but otherwise appear consistent with a spherical distribution. 
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Fig. 10. — Globular cluster luminosity functions for the red and blue samples. A Gaussian profile 
with a turnover magnitude appropriate for the full sample has been superposed and arbitrarily 
scaled in the vertical direction. 
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